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ABSTRACT 

Direction sensitive direct detection of Weakly Interacting Massive Particles (WIMPs) as dark matter would 
provide an unambiguous non-gravitational signature of dark matter (DM). The diurnal variation of DM signal 
due to earth's rotation around its own axis can be a significant signature for galactic WIMPs. Because of 
particular orientation of earth's axis of rotation with respect to WIMP wind direction, the apparent direction 
of WIMP wind as observed at a detector can alter widely over a day. In this work we calculate the directional 
detection rates with their daily and yearly modulations in earth-bound dark matter experiments considering 
detailed features of the geometry and dynamics of the earth-sun system along with the solar motion in galactic 
frame. A separate halo model namely the dark disc model other than the usual standard halo model for dark 
matter halo is also considered and the results for two models are compared. We demonstrate the results for two 
types of gas detectors namely DRIFT (target material CS2) and NEWAGE (target material CF4) that use Time 
Projection Chamber techniques for measuring directionality of the recoil nucleus. The WIMP mass and recoil 
energy dependence of the daily variation of event rates are computed for specific detector and the sensitive 
ranges of mass and recoil energies for the considered detector are probed. 
Subject headings: Dark Matter detection, Diurnal and Annual modulations 



1. INTRODUCTION 

Dark matter is a non-luminous, non-baryonic, pervasive 
fluid extending in and around the bound astronomical 
systems like galaxy and galaxy clusters. This hypothetical 
proposition emerged from an attempt to reconcile astro- 
nomical observations with Newtonian dynamics. Indirect 
indications of existence of dark matter is perceived through 
their putative gravitational interactions in several cosmo- 
logical observations like rotation curves of spiral galaxies, 
the gravitational micro-lensing, observations on Virgo and 
Coma clusters (McLaughlin 1998; Lokas ef al. 2003) bullet 
clusters (Bradac 2009) etc. These observations indicate the 
existence of roughly spherical halo of dark matter around the 
visible galaxies. Most of the constituents of dark matter is 
believed to be non-relativistic (cold dark matter) with Weakly 
Interacting Massive Particles (WIMPs) as the plausible can- 
didates. The direct detection of the dark matter in terrestrial 
detectors involves the scattering of the target nuclei in a 
detector when dark matter interacts with them. The resulting 
recoil nuclei deposit their energy in the detector to produce 
dark matter signal in the detector. Dedicated detectors are 
operational in search of any non-gravitational interaction of 
dark matter with ordinary matter to provide more insight into 
the nature of dark matter and its interactions. The dynamics 
of rotation of the Milky Way galactic disc through the halo 
of dark matter causes the earth to experience a wind of 
WIMPs apparently flowing along a direction opposite to 
that of the motion of the solar system relative to the dark 
matter halo. There are several ongoing/proposed exper- 
imental efforts (Ahmed et al. 2009; Bernabei et al. 2004; 
Bernabei et al. 2008; Angle et al. 2008; Aprile et al. 2009; 
Angloher et al. 2008; Angloher et al. 2005; 
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Cozzini et al. 2005; Ninkovic et al. 2005; 

Burgos et al. 2007; Snowden et al. 2000; Ohnuki et al. 2001; 
Broniatowski et al. 2009; Rozov et al. 2010; 

Sanglard et al. 2005; Nishimura et al. 2008; 

Tanimori et al. 2004; Rubbia 2006; Beltran 2008; 
Billard et al. 2009; Richer et al. 2009; Ahlen et al. 2010; 
Moulin et al. 2005; Cline et al. 2003; Cline et al. 2001; 
Bovine? a/. 1999; Cebrian et al. 1999; Cebrian et al. 2000; 
Cebrian et al. 2001; Battat et al. 2009; Sciolla et al. 2008) 
to directly detect these galactic WIMPs by measuring the 
energy deposition by the WIMPs through their scattering 
off the detector nuclei. Unambiguous identification of a 
WIMP signal is severely challenged by presence of similar 
background signatures. 

An annual modulation is expected (Drukier et al. 1986; 
FreeseefaZ. 1988) in the detection rates of WIMPs in the 
direct detection experiments. This seasonal variation is an 
effect of earth's revolution around the sun in course of which 
the component of earth's orbital velocity parallel to the WIMP 
flow direction annually varies causing the earth to encounter 
different WIMP fluxes at different times of the year. But the 
relative motion between the sun and the earth is much slower 
compared to the motion of solar system in galactic halo. 
Owing to this, the annual modulation is expected to be only 
a few percent of the WIMP detection rate and it is extremely 
difficult to extricate it as a positive signature of dark matter 
from seasonal variation of background rates. On the contrary, 
observation of directional anisotropy of the WIMP wind 
on earth is a signature which can hardly be mimicked by 
any other backgrounds and is potentially more powerful in 
providing unambiguous signature of galactic WIMPs. There 
are dedicated experiments with gas detectors with a Time 
Projection Chamber (TPC) (DRIFT (Burgos et al. 2007; 
Snowden et al. 2000; Ohnuki et al. 2001), NEWAGE 
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(Nishimura et al. 2008; Tanimori et al. 2004), MI- 
MAC (Billard et al. 2009; Richer et al. 2009; 
Ahlen et al. 20 1 0; Moulin et al. 2005), DMTPC 
(Battat et al. 2009; Sciolla et al. 2008)) or organic crys- 
tals (Sekiya et al. 2003) as detectors, aiming to measure 
the directional (angular) distribution of WIMP induced 
nuclear recoils. The technique involved in the directional 
measurements is the 3-D or 2-D reconstruction of nuclear 
recoil tracks with a good spatial resolution followed by 
the determination of sense of the recoil direction (dis- 
crimination between the head and tail of the recoil track). 
The directional distribution of the nuclear recoils which 
carries the imprints of incoming WIMP directions can be 
exploited to disentangle WIMP signals from backgrounds. 
There is a plethora of papers where comprehensive studies 
of various aspects of directional detection of dark matter 
has been performed (Green et al. 2008; Green et al. 2007; 
Copi et al. 2001; Lisanti et al. 2009; Cerdeno et al. 2010; 
Vergados et al. 2007; Gondolo 2002; Buckland et al. 1994; 
Lewinefa/. 1996; Alenazi et al 2008; Copi et al. 1999; 
Copi et al. 2003; Brhlik et al. 1999; Gelmini et al. 2001; 
Belli et al. 2002; Fornengo et al. 2003; Green 2003; 
Morgan et al. 2005; Billard et al. 2010; Vergados et al. 2010; 
Grignon et al. 2010; Creswick et al. 2010). However, be- 
cause of earth's rotation about its axis which has a particular 
orientation relative to the direction of incoming WIMPs, 
the average direction of the WIMP wind with respect to 
an observer fixed on earth, changes over a sidereal day. 
Such modulations in the apparent WIMP direction leads to 
daily fluctuations in the directional (angular) distribution 
of WIMP signal observed at a static detector on earth. 
Observation of such diurnal variations in the directional 
distribution of the WIMP signal would be a robust signa- 
ture for WIMPs (Spergel 1988; Collar et al. 1992). The 
directional dependence, however, is sensitive to the velocity 
distribution of WIMPs in the halo of dark matter. Generally 
the dark matter halo is considered to be spherical and 
isothermal (Binney et al. 2008; Lynden-Bell et al. 1967) 
with the velocity distribution of dark matter WIMPs with 
respect to the galactic rest frame being Maxwellian which 
is a steady state solution of the collisionless Boltzmann 
equation (Chaudhury et al. 2010). Halo model dependence 
of the directional distribution has been extensively discussed 
in (Alenazi et al. 2008). 

The directional sensitivity is usually realised and presented 
in terms of galactic coordinates where the description is 
independent of daily time. In order to comprehend the 
manifestation of directional sensitivity in the WIMP signal 
as registered in the earth-bound detectors, it is essential 
to describe the directional dependence of recoil rates with 
respect to laboratory (detector) fixed coordinate system. The 
transformation from galactic frame to laboratory frame natu- 
rally invokes temporal variation in the directional description 
of WIMP rates over a sidereal day. In this work we explicitly 
calculate the transformation between galactic and laboratory 
reference frames taking into account details of the earth-sun 
geometry and dynamics and compute the time dependence 
of the directional recoil rates in direction sensitive terrestrial 
dark matter experiments. In order to demonstrate the results 
of the directional rates calculations in the laboratory frame 

1 For a detailed description of directional detection of dark matter see 
(Sciolla 2009). 



we assume the standard halo model of dark matter with a 
Maxwellian velocity distribution for the WIMPs. 

In a previous work (Vergados et al. 2009) also the 
authors have studied the diurnal variation of the direct 
detection of dark matter in directional experiments. In 
(Vergados et al. 2009), they express the WIMP velocity in a 
coordinate system where the polar axis (Z) is in the direction 
of recoiling nucleus. In the present work, however, the 
detector fixed coordinate systems are used and the directional 
detection rates are explicitly computed for two specific detec- 
tors. In this case the direction is specified by the direction of 
recoil nucleus with respect to this detector fixed system. For 
the study of diurnal variation in (Vergados et al. 2009) the 
orientation of a detector is specified in the local frame which 
is given by a specific point in the sky designated by the right 
ascension a and declination <5. On the contrary, in the present 
work we explicitly show, for specific detectors, the actual 
variation of the WIMP detection rates over a day for different 
recoil directions (of the target nucleus) with respect to the 
detector fixed frame of reference. 

It is also worth investigating the impact if other halo model 
different from the usual standard halo model is used for 
these calculations. On the other hand, such investigations 
may be useful to probe the structure of dark matter halo 
from the results obtained in actual experiments. In this 
paper we consider another halo model namely dark disc 
model (Read et al. 2008) and compare the results with those 
obtained for standard galactic halo model. 

The paper is organised as follows. In Sec. 2 we present the 
expressions required to calculate the directional recoil rates. 
Sec. 3 describes the analytic calculations for obtaining the re- 
sults in laboratory frame of reference. In Sec. 4 we compute 
the diurnal and annual variation of directional detection rates 
for two terrestrial dark matter experiments. The response of 
the detector to the daily variation of event rates for different 
ranges of WIMP masses and recoil energies is computed in 
this section. Finally in Sec. 5 we give some concluding re- 
marks. 

2. ANGULAR DEPENDENCE OF THE DIRECT DETECTION RATES 

The differential rate for WIMP induced nuclear recoils per 
unit detector mass is given by (Alenazi et al. 2008) 

dE tfliftn J v>Vn v 

where, the summation index n represents a particular nu- 
clear species of the target nucleus of the detector (e.g. n = 
C,S for CS 2 target in DRIFT detector and n = C,F for CF 4 
target in NEWAGE detector). C„ gives the mass fraction 
of nuclear species n in the target material. In the above 
<r„(E) = E miix {d<j„/dE), (d(7 n /dE) being the differential cross 
section of WIMP-nucleus(n) scattering and E m . dx is the max- 
imum possible energy transfered between a WIMP and a 
nucleus of species n, which can be expressed in terms of 
WIMP speed (v), nuclear mass (M„) and WIMP mass (m) as 
Easa = 2filv 2 /M„, with [i n = mM n /(m+M„). In this work the 
local dark matter density p is taken to be 0.3 GeV cm" 3 and 
the detection efficiency (£ (E) for recoil energy E) to be 100%. 
Following (Alenazi et al. 2008) (and references therein), the 
Radon transformation of the velocity distribution /(v) enables 
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one to write the differential rate for WIMP induced nuclear 
recoils per unit solid angle per unit detector mass as 

, P , dR a ,, = Y] -r-^-C n f n ( Vn ,Y/)a n (E)S(E) (2) 
dE dcosOdcp 4-nfipn 

The recoil momentum spectrum for nucleus n is denoted by 
/n(v n ,w) where w is a unit vector along the direction of nu- 
clear recoil and v„ (= c-J 'M„E /2/i|) represents the minimum 
velocity of WIMP required to transfer an amount of energy 
E to the nucleus n (c is the speed of light). The function 
/n(v«,w) is the 3D Radon transform of the velocity distri- 
bution /(v) and is given by /„(v„,w) = J S(y ■ w- v„)f(y)d 3 v 
(Alenazi et al. 2008). For a Maxwellian velocity distribution 
with velocity dispersion <r v for the WIMPs, the recoil momen- 
tum spectrum in the laboratory rest frame is given by 

. " (v„-w-V) 2 

fn(v„ , W) = — r-yr exp 



(27ra2)V2 



2^ 



(3) 



where, V is the average velocity of the WIMPs with respect 
to the detector. The direction of nuclear recoil in three dimen- 
sion is described in a chosen frame of reference by two angles 
9 and tj> as shown in Fig. 1. The recoil energy-integrated di- 




FlG. 1. — Specification of direction of nuclear recoil in WIMP-nucleus 
scattering. 



rectional differential rate is then given by 

dR f ( dR 

dE 

£ih \dEdcos9d(j) 



dcosO dc, 

E± being the detector threshold energy. 



(4) 



3. CALCULATION OF W ■ V 

In this section we present the analytic calculations of the 
direction sensitive term w • V (Eq. 3) in terms of the angles 8 
and 4> that specify the recoil direction w in a terrestrial dark 
matter detector. We describe the direction sensitivity in a lab- 
oratory fixed frame of reference ^Lab, denoted in terms of a 
right-handed system (eLabi, CLab2, CLab3X which are unit vec- 
tors pointing respectively towards east, north and vertical di- 
rection in the sky at the detector location. As shown in Fig. 
2, a particular direction of nuclear recoil, denoted by the unit 
vector w, can be expressed in terms of its components in the 
reference frame S^b as 

w = sin 9 cos <f> eLabi + sin 9 sin (f> eLab2 + cos 9 eLab3 (5) 

The average velocity V of WIMP with respect to detector can 
be written as (Alenazi et al. 2008) 



V = Vwg _ Vsg - Ves - Vle 



(6) 



where, Vwg is the velocity of WIMP with respect to the galac- 
tic centre. In standard halo model Vwg = 0. Vsg is the veloc- 
ity of sun relative to the galactic centre which points towards 
the Cygnus constellation. Ves is the velocity of the centre 
of mass of the earth relative to the sun and Vle is the veloc- 
ity of the laboratory (detector) relative to the centre of mass 
of the earth. In the following subsections we present explicit 
calculations for the quantities w • Vle, w ■ Ves an d w • Vsg- 
In obtaining the analytical expressions we have neglected the 
small eccentricity (0.0167) of the Earth's orbit around the sun. 




towards 
east 



FIG. 2. — Illustration of Laboratory fixed reference frame SLab and the 
direction of nuclear recoil. 



The directional dependence of the observed recoil rates 
convolves the WIMP velocity distribution in the dark mat- 
ter halo and average velocity of the WIMPs with respect to 
the detector. The direction sensitivity of recoil rates are de- 
scribed here, in terms of the two angles 9 and <fi, defined in 
the rest frame of the detector. Hence our computation of ex- 
pression in Eq. 4 directly gives the angular dependence of 
rates as actually expected to be observed in a terrestrial de- 
tector. In a detector-fixed frame of reference the velocity of 
the WIMP wind acquires a time-dependent description due to 
earth's revolution around the sun and earth's rotation about its 
own axis. As a result the observed angular dependence of the 
recoil rates, described in the frame of the terrestrial detector, 
would exhibit both annual and diurnal variations. We compute 
such annual and diurnal variations of the angular distribution 
of the event rates in direction sensitive dark matter detectors. 



3.1. Expression for w • Vle 

The velocity of the detector with respect to the centre of the 
mass of the earth is along the direction of eLabi and can be 
written as 

Vle = | Vle | e La bi = uR® cos a e La bi (7) 

R& being the radius of earth and u> (= 27r/24 rad/hour) is the 
angular speed of earth's rotation about its axis and a is the 
latitude of the detector location. Therefore 

w- Vle = w/?© cos a sin # cos ^ (8) 

3.2. Expression for w ■ Ves 

To find w • Ves we bring in three different frames of refer- 
ences: Seq, S\q and S20 which are chosen as follows. The 
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frame S eq , as shown in Fig. 2, has a origin fixed at the centre 
of mass of earth and is described by a right-handed system of 
orthonormal vectors (e e qi,e e q2,e 3 ). The plane spanned by the 
unit vectors e eq i and e eq 2 defines the plane of the equator with 
e eq i pointing vertically upward direction at equator at 12 mid- 
night. e 3 is a unit vector pointing towards direction of angular 
velocity of earth's rotation and therefore defines the earth's 
axis of rotation. The transformation between the frames ^Lab 
and S eq frame is given by 



/eeqi\ /-sinu;f-sina:coscjf cos a cos a; A /eLabi\ 

e eq 2 = coswt -sinasinojf cos a sin cut eLab2 (9) 
V e3 / V cos a sin a / \e L ab3/ 

where, t is the time elapsed after 12 midnight. S\® and S^q 
are two frames of reference with their origins fixed at the sun 
and are shown in Fig. 3. They are described respectively by 
the two right-handed systems (ei , e2, 63) and (e [ , e2 , e 3 ) where, 
ei, e2 are mutually perpendicular unit vectors in the ecliptic 
plane, with ei pointing towards the centre of mass of earth 
at summer solstice from the sun. e3 is a unit normal to the 
ecliptic pointing towards the orbital angular velocity of earth. 
The axis of earth's rotation for all time of the year is parallel 
to ei-e3 plane making an angle S = 23.5° with the e3 and at 
summer solstice its tilt is towards sun. The mutually orthogo- 
nal unit vectors e'j and e2 define the plane of the equator. The 
transformation between frames Siq and S20 is then given by 




FIG. 3. — Illustration of reference frames Siq, 520 and Seq (described in 
text) 




cos 8 - sin 8 

1 
sin 5 cos 6 




(10) 



A unit vector along the radius vector of earth with respect 
to the sun at a time T of the year (counted from the day of 
summer solstice) is given by 



e r = cos fir ei + sin fir e2 
= cos fir cos <5 ej + sinilT ez~ cos sin 8 e 3 



(11) 



where, fi = (2ir radians/year) is the orbital angular speed of 
the earth about sun. Since e eq i is along the vertical upward 
direction at equator at 12 midnight it is parallel to the line 
of intersection of two planes: the plane containing the axis 
of the rotation of earth (e 3 ) and the radius vector joining sun 
to centre of earth and the plane of the equator. The normal 
to the two planes: the plane containing the vectors e 3 & e r 
and the plane of the equator respectively are e 3 x e r and e 3 . 
The perpendicular to both of these normals define the line of 
intersection of the two planes and along which the unit vector 



^e q l : 



(e 3 x e r ) x e 3 



(e 3 x e r ) x e 3 | 
cos fir cos 8 e'j + sin fir e2 

V cos 2 fir cos 2 8 + sin 2 fir 
The unit vector e eq 2 is then obtained as 

Ce q 2 = C3 x e eq i 

-sin e'j + cosSircosJ e2 



(12) 



(13) 



Vcos 2 fir cos 2 5 + sin 2 QT 
Therefore the transformation between reference frames 52© 
and S eq is given by 

/ycosfircos<5 -ysinfir 0\ /e eq i\ 

ysinfir ycosfircosi50 J I e eq 2 I (14) 
ij V e 3 J 




where 



1 



y = 



(15) 



\J cos 2 fir cos 2 S + sin 2 fir 
Using Eqs. (9,10, 14) we can write the transformation equation 
between the frames Siq and ^Lab as 

(eA /cosi50-sin(5\ /ycosilT cosS -ysinfir 0\ 
e 2 = 1 ysinfir ycosfir cos<50 

e 3 / Vsin^O cosS J \ ij 

C- sin ujt - sin a cos tot cos a cos ut \ / eLabi\ 
coswf - sin a sin cut cos a sin cut eLab2 (16) 
cos a sin a / \eLab3/ 

A unit vector along the direction of instantaneous orbital ve- 
locity of earth with respect to sun at time T of the year is given 
by ' 

eES=~sinfir ei + cosfir e2 (17) 

Using the transformation in Eq. 16 we can express it in the 
laboratory frame of reference as 

eES=3 ; ^ cos S cos ujt - -sin2firsin 2 <5sinu;f^ eLabi 
+ ( cos a sin 8 sin fir - y sin a cos 8 sin cot 

y 2 \ 

--sinasin2firsin ^coswfj eLab2 
+ ( sin a sin 8 sin fir + y cos a cos 8 sin tot 

y t r \ 

+-cosasin2fir sin Sco&ojtJ eLab3 (18) 
Finally we get 
w Ves = |Ves| (w-e E s) 



= |V ES ] 



y I cosjcoso;/- 



— sin2f2rsin 2 <5sinajf i \inf7cos. 



+ ^ cos a sin 5 sin Q.T -y sin a cos 8 sin ujt 

y 2 \ 

— sinasin2Slrsin Scostvt sin sin 
2 / 

+ ^ sin a sin 8 sin QT + y cos a cos 8 sin cot 
+— cosasin2f2r sin 2 8 cos cut ) cost 



(19) 



where, we have taken the magnitude |Ves| of earth's orbital 
speed as 30 km/s (Plank 2008) 
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FIG. 4. — Left panel: Relative orientation of ecliptic and plane of equator 
at the position of vernal (March) equinox. Right panel: Declination and right 
ascension of Cygnus direction or direction of velocity of sun with respect to 
galactic centre has been illustrated. 

3.3. Expression for w Vsg 

The velocity Vsg of the sun with respect to the galactic 
centre happens to point towards the direction of Cygnus con- 
stellation. The orientation of Cygnus constellation in celestial 
coordinate system is given by a declination of 9 C = 42° and 
a right ascension of 20.62 hours (Ridpath et al. 2007). This 
implies that the vector Vsg makes an angle 90° - 9 C with the 
axis of rotation of earth (e' 3 ) and its projection on the plane of 
the equator subtends an angle <fi c = 20.62/i ~ 309° when mea- 
sured (anticlockwise) from the direction in sky where the sun 
crosses the celestial equator at vernal equinox. This direction 
is defined by e2- Therefore a unit vector along the direction of 
Vsg can be written in the frame S20 : (ej, ^2,^3) as 

(20) 



„_ .. „ _ AVJ . v _ 1; ~„ 

esc = _ cos 9 C sin <fi c e[ + cos 9 C cos <f> c e2 + sin 9 c e' 3 



The transformation matrix between the frames of reference 
S20 : (ei,e 2) e3) and S Lab : (e Labl ,e Lab2 ,e Lab3 ) can be obtained 
using Eqs. (9,14) as 

/ e i\ /ycosJ7rcos<5 -ysinttT Q\ 
ysinflT ycosflT cosSO J 
1/ 
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(— sin tot - sin a cos tut cos a cos wA / eLabi \ 
coswf -sin a sin ut cos a sin ut eLab2 
cos a sin a / \e L ab3/ 

(21) 

Exploiting the above equations, esc can be expressed in terms 
of the laboratory frame of reference as 

esc = (~A sin ut + B cos ut) eLabi 

- sin a(A cos uit + B sin uiT) + cos a sin 9 C eLab2 

+ (cosa(Acoswf +fisincjr) + sinasin6> c eLab3(22) 

where, 

A = y cos 9 C (- sin <fi c cos ttT cos 6 + cos (f> c sin ftT) (23) 
B = y cos <9 c (cos cp c cos UT cos 8 + sin <f> c sin QJ) (24) 

Thus we have 

w Vsg = I Vsg I (w-e S c) 

= I Vsg 1 1 sin 9 cos cf>{— A sin ut + B cos ut) 

+ sin # sin </> -sinct(Acosa;f+Bsina;r) + cosasin6' f 

+ cos0 (cosa(Acoswf +Bsinwr) + sinasin6' c M.25) 



Although the physical range of |Vsg| is given by 170 km/sec 
< |V SG | < 270 km/sec (90 % C.L.) (Leonard et al. 1990; 
Kochanek 1996; Bovy et al. 2009), in the present work we 
consider the central value of |Vsg| as 220 km/s. 

4. DIURNAL AND ANNUAL VARIATIONS OF DIRECTIONAL 
DETECTION RATES 

In this section we calculate the diurnal and annual varia- 
tions of the direct detection rates with respect to the observed 
direction of recoil nuclei with the formalism described in the 
previous sections. In order to estimate such detection rates, 
the detectors are considered to be direction sensitive, i.e. the 
directionality of reconstructed tracks of the recoil nuclei are 
apparent from the observed signals (i.e. the recoil tracks will 
have no head-tail ambiguity). We also assume the detectors to 
be 100% efficient in the present work. From Eq. 2 it is evident 
that, these rates depend on the target materials and location 
(latitude) of the detector. For order of magnitude estimations 
of such rates therefore, the only information of detector prop- 
erties in this calculation are given in terms of the mass and 
atomic numbers of the target nuclei and the mass fraction (C„) 
of each type of target nucleus (in case of diatomic or mul- 
tiatomic targets), as required to compute Eq. 2. In this work, 
we consider two directional detectors namely DRIFT (target 
material CS2) (Burgos et al. 2007; Snowden et al. 2000; 
Ohnukief al. 2001) and NEWAGE (target material CF 4 ) 
(Nishimura et al. 2008; Tanimori et al. 2004). 

The DRIFT detector at Boulby underground laboratory 
(54°33'N,0°49'W) is designed to detect WIMPs via nuclear 
recoils within the fiducial volume of the detector. The detec- 
tor has the ability of not only counting nuclear recoils but also 
determining the direction of nuclear recoils by reconstructing 
three dimensional tracks of nuclear recoils. The DRIFT II 
detector is a Time Projection Chamber (TPC) containing 
gaseous CS2 at a pressure of 40 torr. It consists of two back 
to back TPCs on either side of a central photocathode. The 
ionization electrons resulting from the scattering of WIMPs 
off the C or S atoms of CS2 gas get captured by the elec- 
tronegative CS2 molecules to form CSj ions. These ions drift 
in the electric field towards the anode plate where a multiwire 
proportional chamber consisting of a two dimensional grid 
of wires are used as readouts. The CS2 ions experience 
only thermal diffusion allowing a 3D reconstruction of the 
recoil track range and direction. The NEWAGE is also a 
direction sensitive dark matter search experiment at Kamioka 
underground laboratory (36°25'N, 137°18'E). It uses a three 
dimensional gaseous tracking detector (Micro-TPC) which 
has Micro Pixel Chamber as its readout. The gaseous micro- 
time projection chamber is filled with CF4 at 152 torr and is 
read by a two dimensional fine pitch imaging device called 
/i— PIC consisting of orthogonally formed cathode and anode 
strips. The details of detector fabrication has been described 
in (Nishimura et al. 2009). These kind of techniques are used 
to reconstruct the WIMP induced nuclear tracks in both two 
and three-dimensions which then enable one to determine the 
directionality of the recoil track. The recoil energy threshold 
for the DRIFT experiment can be ~ tens of keV (Plank 2008). 
The NEWAGE experimental collaboration has mentioned a 
recoil energy range ~ 100-400 keV for their WIMP-search 
analysis. The lower limit of this is set by analysing the 
background gamma rejection energy (Miuchi et al. 2007). 
Since the WIMP induced events will be small it is required to 
have very low threshold for the dark matter detectors albeit 
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it increases the unwanted backgroud and unwanted neutron 
signals. The expected WIMP induced recoil energy range 
is ~ 1 - 50 keV. Shielding etc. are required with specific 
techniques to reduce the unwanted signals/backgrounds. In 
this work we have presented the results with zero threshold 
energy for both DRIFT and NEWAGE. 

Owing to the apparent directionality of WIMP wind an 
anisotropy is expected in distribution of nuclear recoils. 
The directional rates are also expected to exhibit a daily 
modulation due to rotation of the earth. Such signals of daily 
modulation can provide an effective background discrimi- 
nation and thereby producing powerful signature of galactic 
dark matter WIMPs. 

There are two dominant types of WIMP-nucleus elastic 
scattering - (i) spin-independent interaction, where WIMP 
interacts with the nuclear mass and (ii) spin-dependent 
interaction, where the WIMP interaction with the nuclear 
spin in considered. The WIMP-nucleus scattering cross 
section can have contributions from either or both of the two 
interactions. In case of DRIFT detector with target material 
CS2, only spin-independent interaction contributes to the 
total scattering cross sections as both C and S nuclei have 
zero spin at the ground state. On the other hand, since F 
nucleus has a ground state spin 1 /2, the NEWAGE detector 
with target material CF4 is sensitive to both the interactions. 
But the contribution from the spin-dependent part comes 
out to be ~ 2-3 orders of magnitude higher than that from 
spin-independent part. Henceforth for NEWAGE we show 
the event rates with spin-dependent cross sections only. 

For a detector with target material CS2, the spin- 
independent directional recoil rate (as given in Eq. 4) can be 
written as (Alenazi et al. 2008) 
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In the above, a is the WIMP-proton cross section, n p is the 
WIMP-proton reduced mass in GeV and f^(E,v/) is given 
by 

f^(E , w) = C C A 2 C /c(v c , w) T c + CsA 2 Js(v s , w) F s • (27) 

The function /c,s(vc,SjW) can be evaluated using Eq. 3. Fc, 
J-"s are squares of the nuclear form factors (F(E)) given 
by F(E) = exp(£/2£'o) where E is the energy transferred 
from WIMP to the target nucleus, Eq = 3/ (Im^R 2 ^), m# be- 
ing the nuclear mass and Rq being the radius of the nu- 
cleus (Jungman et al. 1996). The form factor is more ac- 
curately given as (Jungman et al. 1996; Engel 1991) F(E) = 
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), where the momentum transferred q ■■ 
5s 2 ) 1 / 2 , R~ l.2A l IHm, where j x 



ii = \k~ — js~)'' ", n ~ i.z/i -'-mi, wnere j\ is a 
spherical Bessel function, s = lfm and A is the mass number 
of target nucleus. The form factor depends on the distribution 
of mass and spin within the nucleus and also on the type of 
the WIMP-nucleus interaction. For very small energy trans- 
fer E, relevant for our calculation, the form factor ~ 1 . Ac and 
As are the mass number of carbon and sulphur nuclei respec- 
tively. Cs and Cc are the mass fractions of carbon and sulphur 
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FIG. 5. — These plots shows distribution of energy integrated WIMP- 
nucleon spin-independent event rates in CS9 target (DRIFT) detector over 
8 and <j> at different time of the year and at different instants of the day 
(described in the text) . Color contour shows event rates in units of 
[events/kg/day s/sr] . Plots are done with a cross section of lfr 44 cm 2 and 
a dark matter mass of 60 GeV. 



respectively and are given by 
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(28) 



In the above M c (~ 11.26 GeV) and M s (« 30.04 GeV) are 
masses of carbon and sulphur nuclei respectively. 

The spin-dependent direction recoil rate for a detector with 
target material CF4 is given by the same expression as Eq. 26 
but with /jjj replaced by /JJ? which is given by 
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where CV is the mass fraction of the fluorine nucleus given by 

4M f 



C F = - 



4M F +M C 



(30) 



Mf being the mass (w 17.84 GeV) of the fluorine nucleus. 
The factor (0.647) comes from spin of fluorine nucleus and 
has been discussed in (Alenazi et al. 2008). 

The angular dependence of the spin-independent rates 
(number of events per day) per kg of target material CS2 of 
DRIFT detector are calculated using Eqs. (2-28) . Similar 
calculations for spin-dependent rates with NEWAGE detector 
(CF4) are performed using Eqs. (2-26, 29-30). The results for 
the spin-independent case (DRIFT) are shown in Fig. 5. The 
directional rates plotted here are calculated with a benchmark 
value of WIMP-nucleon cross section of 10~ 44 cm 2 and a 
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FIG. 6. — Same as Fig. 5 but for spin-dependent events in CF4 target 
(NEWAGE) 

WIMP mass of 60 GeV. Each plot shows the variation of 
event rates with 9 and cf> as defined in Sec. 3. They are shown 
in Fig. 5 using a colour code symbol. A colour corresponding 
to a point in the 6 — <p plane for any single plot in Fig. 5, 
therefore indicates the detection rate for the directionality 
given by the 6 — <f> coordinates of that point. The colour 
codes are chosen such that the increase of detection rates 
can be represented by the change of colour from green to 
white in the order green - deep blue - light blue - white. The 
annual and diurnal variations of the directional dependence 
of the event rates observed in the detector are illustrated 
by presenting the plots for different values of the set (T,t), 
where T and t represent a time of the year and a time in a day 
respectively. We show the results for 16 representative sets 
of (T,t) constructed out of four different times of the year 
namely summer solstice, autumnal equinox, winter solstice 
and vernal equinox and for four different times of the day - 
0.00 hrs (midnight), 6.00 hrs (6 A.M.), 12.00 hrs (noon) and 
18.00 hrs (6 P.M.). The plots of a given column represent the 
annual variation of the directional rates observed at a partic- 
ular time of the day. The plots of a given row, on the other 
hand, represent the diurnal variations of the directional rates 
observed at a particular time of the year. For example, from 
Fig. 5 we see that, for (0,0) = (100° , 230°) the directional 
rate is 5.5 x 10~ 6 events/kg/day s/sr at 0.00 hours on the day 
of summer solstice (left corner plot of Fig. 5). As one goes 
down along the first column of plots (annual variation) in 
Fig. 5 this rate changes to 4.5 x 10" 5 , 1.2 x 10" 4 , 2.9 x 10" 5 
in units of events/kg/days/sr for autumnal equinox, winter 
solstice and vernal equinox respectively at 0.00 hours. The 
representative colour for the rates also changes accordingly. 
This rate also shows the variations in different times of a day 
for the chosen value of 8,<f> as one moves row-wise from 
the left corner plot of Fig. 5. In this particular example the 
rate varies from 5.5 x 10~ 6 events/kg/days/sr to 6.5 x 10~ 5 , 
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FIG. 7. — upper panel: ^ (spin-independent) vs <j> at different values of 9 
marked at the inset of 2 nd column. Lower panel: ^ (spin-independent) vs 9 
at different values of <fi marked at the inset of 3 rrf column. We have chosen a 
WIMP-nucleon cross section of 10 -44 cm 2 , dark matter mass of 60 GeV. The 
target material is CS2 (DRIFT detector) 
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FIG. 8. — Same as Fig. 7 but for spin-dependent interactions at a CF4 in 
NEWAGE detector. 



1.5 x 10~ 4 , 2.7 x 10" 5 (in units of events/kg/days/sr) that 
correspond to three different times of the day namely 6.00, 
12.00 and 18.00 hours respectively on the day of summer 
solstice. As is evident from Fig. 5, these daily variations of 
the rates are also indicated by the change of representative 
colours. Similar observations can be made for plots in any 
other chosen rows or columns. In Fig. 6 we present the same 
plots as in Fig. 5 but for spin-dependent events in a CF 4 
target of NEWAGE detector. The pattern of variation of event 
rates with 9 and <\> differs in the two figures (Figs. 5 and 
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6). From Eqs. (27) and (29) we expect amplitude variation 
between the spin-dependent and spin-independent case. 
Again from Eq. (25) along with Eqs. (2) and (3) one sees 
that the directional rates also depend on the latitude (a) of 
the detector location 2 . This is reflected in Figs. 5 and 6. The 
diurnal variations of angular dependence of the directional 
rates are also shown for both the cases of spin-independent 
and spin-dependent WIMP-nucleus interactions in Fig. 7 
and Fig. 8 respectively for more clarity. In these two figures 
we have shown the variations of the directional event rates 
dR/d(cos9)d<fi with 9 at different values of <j> (upper panel) 
and with cf> at different values of 9 (lower panel). The plots 
are shown at four different times of a given day (the summer 
solstice). From Fig. 7 it is evident that, for CS2 detector, 
the range of directional variations (variations in 9, <j> ) of the 
spin-independent rates lies between ~ (5 x 10~ 6 -2 x 1CT 4 ) 
events per day per unit solid angle per kg of detector mass. 
The corresponding range for spin-dependent rates in a CF4 
detector is (2 x 10" 9 -2 x 10" 7 ) events/day/kg/sr ( Fig. 8 ). 
We also show in Fig. 9 the daily variation of 9,<fi- integrated 
rates for DRIFT detector for WIMP mass of 60 GeV. The 
large daily variation is also apparent from this figure. 
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FIG. 9. — Daily variation of 8,(j>- integrated rates for DRIFT detector for 
WIMP mass of 60 GeV 
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FIG. 10. — The number of nuclear recoils in two solid angles with dif- 
ferently chosen directions namely, north and east in the laboratory frame of 
reference as a function of time of a day. The plots are shown at the day of 
summer solstice. 
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FIG. 1 1. — Schematic diagram to show the choice of solid angles as de- 
scribed in the text. 



It is also evident from Figs. 5-8 that calculated daily 
and annual variations of detected number of events within a 
given solid angle is sensitive to the chosen direction about 
which the solid angle is considered. For example, event 
number calculated for DRIFT around the direction specified 
(in the Laboratory frame) by 9 = 90° and = 300° has a 
larger variation over a day and also over the year. This 
can be understood from Fig. 5 as the representative colour 
associated with the point exhibits a span of variations ranging 
from green to white over different times of a day (or a year). 
On the contrary, number of events calculated around the 
direction specified by 9 = 20° and <fi = 25° shows mild daily 
and annual variations. This is realised from Fig. 5 since over 
different times of a day (or a year) the point lies mostly in 
the green region of the 9 - <f> plane. The daily variations of 
direct detection rates for three chosen directions of recoil 
tracks are shown in Fig. 10 for illustrative purposes. In 
Fig. 10 we plot the calculated number of events within a 

2 DRIFT is at a latitude a = 54°33'N, for NEWAGE a = 36°25'N 



solid angle considered about two different chosen directions 
namely north and east as a function of different times of 
a sidereal day. In this case we have chosen this day to be 
the day of summer solstice. The solid angles are so chosen 
that the semi-vertical angle of the cone around the chosen 
direction is 15°. It reveals larger daily variations along east 
direction in comparison to the northern direction for this case. 
Experiments dedicated to directional dark matter detection 
may therefore look for some preferential directions as 
obtained from such theoretical studies for the daily variations 
calculated in the Laboratory frame of reference, along which 
the substantial daily variation is expected. 

The diurnal variations in the directional distribution of 
WIMP signals can be better illustrated by studying the daily 
fluctuations of a dimensionless quantity defined as the ratio 
of directional recoil rate in a direction Wi(6\,<fii) to that in a 
direction W2(#2, 4>2)- In a practical scenario this can be esti- 
mated by considering ratio of number of events for two solid 
angles £l\ and around the two specified directions (wi and 
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Wa) which can be defined as 

It 
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(31) 



Such variations, estimated for the DRIFT detector, are 
shown in Fig. 12 where we have presented the daily variation 
of the quantity r[wi,w 2 ] on the day of summer solstice for a 
representative set: (Wi,Wi) chosen as {east (6 = 90°, 4> = 0°), 
north(6> = 90°, <j) = 90°)}. with the semi vertical angle of the 
cones (that define the solid angles considered here) taken to 
be 15°. The choice of the solid angles are illustrated in Fig. 
1 1 where two choices of solid angles are shown - Our calcu- 
lations show that the ratio of detected events corresponding to 
the eastward recoil direction to that corresponding to the re- 
coils along the north varies in the range ~ 0.001 -21 over a 
sidereal day on the day of summer solstice. 
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FIG. 12. — The ratio of number of nuclear recoils in two solid angles with 
differently chosen directions in the laboratory frame of reference as a function 
of time of a day. The plots are shown on the day of summer solstice. We have 
presented the ratio for a representative set :{w\ , wi): (east, north) 



It may be of interest to investigate the daily variation of 
detection rates as functions of WIMP masses and recoil 
energies. In Fig. 13 we show the energy integrated rates for 
different WIMP masses at different time of the day. These 
rates are plotted for two detectors - DRIFT and NEWAGE - 
considered here. The overall yield for each of the detectors 
is maximum at WIMP masses around 20 GeV. The daily 
variation also dominates for WIMP masses at around 20 GeV 
for both the detectors considered here. We also calculate the 
daily variation of rates with recoil energy for a fixed WIMP 
mass of 60 GeV and the results are shown in Fig. 14. From 
Fig. 14 one sees that the overall yield decreases with the 
increase in recoil energy although the daily variation is more 
prominent for higher recoil energies. We also investigate the 
sensitivity of mass ranges to the daily variations for different 
ranges of recoil energies. For demonstrative purposes we 
show the results for DRIFT detector only. In Fig. 15 the 
results are plotted for five fixed values of WIMP masses 
ranging from 10 GeV - 200 GeV. All the rates shown in this 
figure are integrated rates over 8 and <j) and over the con- 
sidered energy ranges. The daily variations are investigated 
for five recoil energy ranges: 0-10 keV, 5-15 keV, 10-20 
keV, 15-25 keV, 20-30 keV. Beyond 30 keV the rate falls 
down to very low value. It is evident from Fig. 15 that the 
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FIG. 13. — Daily variation of detection rates as a function of WIMP mass 
for DRIFT and NEWAGE detectors. 
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FIG. 14. — Daily variation of detection rates as afunction of recoil energy 
for DRIFT and NEWAGE detectors for a given WIMP mass of 60 GeV. 
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FIG. 15. — Daily variation of detection rates as afunction of recoil energy 
and WIMP mass for DRIFT detector. 



overall rate dominates for lower recoil energy range (0-10 
keV) all along the WIMP mass range considered. The nature 
of daily variations of the detection rates are also similar for 
the whole mass range in the recoil energy band of - 10 
keV. It is also observed from Fig. 15 that for higher recoil 
energy ranges, the yield for lower WIMP masses gradually 
decreases but the magnitude of daily variation increases. It 
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FIG. 16. — Comparison of diurnal variation of directional detection rates 
for dark disc model and standard halo model for dark matter for two sets of 
(6, <j>) on the day of summer solstice. 
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FIG. 17. — Diurnal variation of directional detection rates for dark disc 
model with different dark matter mass and recoil energy range. The plots are 
shown for a given 0(45°), 0(120°) and on the day of summer solstice 



appears from the present calculations that for the particular 
detector considered, lower recoil energy range is more suited 
for investigating dark matter if the WIMP masses are having 
a wide range of values (from very low to very high). But if 
WIMP masses are limited to only lower values then the daily 
variation will be more prominent for higher recoil energy 
ranges. 

Diurnal variations of detection rates in alternative galac- 
tic halo models: We also calculate the diurnal variations of 
the directional detection rates considering a different galac- 
tic halo model. Such consideration enables us to compare 
the results obtained from the standard halo model with that 
from this new model. The actual experimental results when 
compared with such calculations may throw light on the na- 
ture of dark matter halo. To this end, we consider a model 
proposed by (Read et al. 2008) which is a combination of the 
usual standard halo model and a proposed dark disc that coro- 
tates with the galactic disc. In the dark disc (DD) model it is 
assumed that the DD kinematics matches with that of Milky 
Way's thick disc (Bruch et al. 2008). This in turn, gives a ro- 
tation lag represented by Vi ag with respect to Vsg described 
earlier in Eq. 6 along with a different velocity dispersion. A 
Maxwellian velocity distribution for dark matter is also con- 
sidered for dark disc model (Bruch et al. 2008) with Vsg re- 
placed by Via e and in this calculation a dispersion of 50km/s 
is adopted following (Bruch et al. 2008). In Fig. 16 we 
show the daily variation of rate on a particular day (summer 



solstice) for dark disc model and compared them with that of 
the standard halo model (SHM) for two sets of values of 8, <fi. 
The rates from the combined model of DD and SHM are also 
displayed in Fig. 16. From Fig. 16 one sees that the yield in 
case of standard halo model is more than that for dark disc for 
the chosen value of dark matter mass. The difference in the 
yield for the two cases varies when observed in different 6,4> 
directions. The magnitude of daily variations also appears 
to differ marginally for the two models. However we have 
checked that this difference between the two models varies 
for different recoil energies and dark matter mass. In Fig. 17 
the diurnal variation for two chosen dark matter masses are 
shown for different recoil energy ranges considering only the 
dark disc model. One of the chosen dark matter mass is in the 
high mass range (200 GeV) and the other one is chosen at 60 
GeV as before. From Fig. 17 one sees that if only dark disc 
model is considered, the yield is much higher for the higher 
dark matter mass and low recoil energy range implying that 
the manifestation of the dark disc effect is more prominant 
for low recoil energy and large WIMP mass. 

5. CONCLUSION 

Directional detection of dark matter in direct detection 
experiments not only provides unambiguous signature of 
galactic WIMPs but also serves as a potentially powerful 
probe to the structure and dynamics of dark matter halo. 
The directional detection techniques involve reconstruction 
and head-tail discrimination of tracks of nuclear recoil, 
induced by the scattering of WIMPs off the detector nuclei. 
Directional description of the WIMP induced nuclear recoil 
rates in terms of angular coordinates defined and realised in 
the laboratory frame of a detector exhibits temporal variation 
over a sidereal day due to diurnal motion of the earth. It 
would also have an annual modulation owing to yearly 
variations of the component of earth's orbital velocity along 
incoming WIMP direction. Analytical details of the annual 
and diurnal variations are determined by the dynamics and 
geometry of the sun-earth system and the location (latitude) 
of the detector on the earth. Proper transformation of coor- 
dinate systems are also needed to obtain the results in terms 
of laboratory frame of reference. This is very crucial for any 
detector based prediction of such directional modulations of 
dark matter signals. The directional rates of WIMP signal 
events in directional detectors sensitive to spin independent 
WIMP-nucleus interactions (CS2 target) are estimated for all 
possible recoil directions and they are found to lie within the 
range ~ 10~ 6 - 10~ 4 per day per steradian per kg of target 
mass of the detector. A WIMP mass of 60 GeV and WIMP- 
nucleon cross section of 10 _44 cm 2 have been assumed for all 
the calculations. Similar estimations for directional detectors 
with dominant sensitivity to spin-dependent WIMP-nucleus 
interactions gives the corresponding range in same units as 
~ 10~ 9 - 10~ 7 . The diurnal variations in the laboratory frame 
description of directional anisotropy of the WIMP induced 
nuclear recoils amount to a wide daily fluctuations (ranging 
from 0.3 to 20) of the ratio of WIMP rates corresponding to 
two chosen orthogonal directions of nuclear recoils (east and 
north). 

Energy integrated rates for different WIMP masses ranging 
from 0-100 GeV show that the daily variation is most 
prominent for WIMP masses around 20 GeV. The study (with 
WIMP mass of 60 GeV) reveals that although the range of 
daily modulation increases with the increase in recoil energies 
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but the drop in overall yield for higher recoil energies makes 
it difficult for a detector to register this modulation. The effect 
of varied recoil energy ranges on the integrated (over 6,<f> 
and over the considered energy range) yield and their daily 
modulations at DRIFT detector for a wide range of WIMP 
masses (10-200 GeV) are also investigated. The results show 
similar daily modulations for the mass range considered for 
lower recoil energy range (0-10 GeV). The overall yield 
is also higher in this recoil energy range. Higher recoil 
energy ranges, although results in the increase in magnitude 
of modulation for lower WIMP masses but the overall yield 
drops down. Hence, if the detector is efficient and precise 
enough to detect the yield corresponding to higher recoil 
energy ranges then for lower values of WIMP masses this 



detector will observe considerable daily modulation of the 
signal (as shown in Fig. 15). 

The results of the analysis are found to differ mostly in 
terms of the total yield when another model for dark matter 
halo namely the dark disc (where one considers a corotating 
disc of dark matter with the galactic disc) model or DD model 
is considered. However, calculations using a combined model 
of DD and the usual SHM shows marginal increase in overall 
yield and also a variation of nature of daily modulations for 
certain recoil directions. 

Acknowledgements: D.M. thanks J D Vergados for some 
discussions. 
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